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We present a dynamic density functional theory (dDFT) which takes into account the advection of the particles 
by a flowing solvent. For potential flows we can use the same closure as in the absence of solvent flow. The 
structure of the resulting advected dDFT suggests that it could be used for non-potential flows as well. We apply 
this dDFT to Brownian particles (e.g., polymer coils) in a solvent flowing around a spherical obstacle (e.g., 
a colloid) and compare the results with direct simulations of the underlying Brownian dynamics. Although 
numerical limitations do not allow for an accurate quantitative check of the advected dDFT both show the same 
qualitative features. In contrast to previous works which neglected the deformation of the flow by the obstacle, 
we find that the bow-wave in the density distribution of particles in front of the obstacle as well as the wake 
behind it are reduced dramatically. As a consequence the friction force exerted by the (polymer) particles on the 
colloid can be reduced drastically. 

PACS numbers: 



I. INTRODUCTION 

The generalization of classical density functional theory 
(DFT) to non-equilibrium states has become a valuable tool 
to study the dynamics of interacting Brownian particles. This 
dynamic density functional theory (dDFT) for the ensemble 
averaged density was proposed recently in 11101 ■ Although 
hydrodynamics is known to play a crucial role in the dynam- 
ics of suspensions, dDFT has also been used to describe col- 
loidal suspensions or polymer solutions, in particular to in- 
vestigate the distribution of solute particles around a strongly 
repulsive potential moving through the solution of particles 
yD. The intention was to model a colloidal particle moving 
through a polymer solution. A similar model has been used to 
study the depletion interaction between two colloidal particles 
moving through a polymer solution 10,111]. In all these stud- 
ies the hydrodynamic flow of the solvent around the colloid 
was neglected and the solvent effectively passed through the 
colloid. A real colloid would displace the solvent as it moves, 
as shown for the case of a small and a large colloid in Fig.[T] 
For a spherical colloidal particle of radius R dragged through 
an unbounded incompressible viscous Newtonian solvent with 
velocity c at low Reynolds number, the flow field u(r) (in a 
frame of reference comoving with the colloid) is given by the 
solution of the Stokes equation i^, 

(1) 

For large distances from the colloid, r ;» i?, the flow field is 
well approximated by u(r) = — c. For solute particles (e.g., 
polymers or other colloidsVwhich only feel this far field, the 
model presented in llH 0, Ht] is a reasonable approximation. 
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This is the case for large solute particles with a radius R. 
Their centers can approach the dragged colloidal particle only 
up to a distance D ~ R + d d, see Fig.flla), and thus will 
feel a flow field u(r) « — c, if one neglects the additional ef- 
fect of the solute particles on the solvent flow, in other words, 
the hydrodynamic interaction between the solute particle and 
the colloid is neglected. This is a reasonable approximation 
for polymer coils but certainly a bad one for solid solute parti- 
cles. But small solute particles of radius d ^ R can get much 
closer to the colloid and feel the distortion of the solvent ve- 
locity field as illustrated in Fig.fTlb). The solute particles will 
be deviated from the colloid by the flow field and this will re- 
duce significantly the bow-wave effect in front of the colloid 
presented in HI and the strenght of the non-equilibrium deple- 
tion force discussed in |0,|5t]. Extremely small solute particles 
would not show this effect at all since they would behave like 
solvent molecules. In this limit, however, the basis of the the- 
ory discussed here, i.e, the description of the solute particles 
as overdamped Brownian particles, is no longer valid because 
it is based on a separation of the length and time scales asso- 
ciated with the solvent molecules and the solute particles. 

In this paper we present a generalization of the dDFT de- 
rived in 111, [jl to the case of Brownian solvent particles ad- 
vected by a flow, thereby incorporating some aspects of the 
hydrodynamics of the solvent into the theory. However, we 
do not model hydrodynamic interactions between the solute 
particles as well as the back-reaction of the solute particles on 
the flow field, e.g., by a concentration dependent viscosity, or 
by a reduced mobility of the solute particles in the vicinity of 
the colloid. However, the latter can be included in a straight- 
forward manner as we discuss in the conclusions in Sec. |IV] 
In the following Sec.|ll]we derive the advected dynamic den- 
sity functional theory using the method described in |0, [Ml- 
In Sec.|lII]we discuss two sample cases, namely ideal solute 
particles and Gaussian solute particles which stress the impor- 
tance of taking into account the solvent flow. 
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FIG. 1: Cross section of the flow field given by Eq. l[T} in a plane par- 
allel to the direction of motion around (a) a small spherical colloid 
and (b) a big one (full circles, radius R). The dashed circles of radius 
D mark the points of closest approach of solute's centers (point in 
the center of open circles) to the colloid. The solute's diameter in its 
mutual interaction is a and, for non-additive mixtures, not necessar- 
ily equal to its diameter 2d = 2 {D — R) (indicated by the dotted 
circle) in the interaction with the colloid. The component of the flow 
field normal to the dashed circle is larger for the small colloid (a) 
than for the large colloid (b). a and d are the same in both figures. 



II. ADVECTED dDFT 

We start with the Langevin equation of an ensemble of N 
advected interacting Brownian particles confined to a finite 
volume V in the overdamped limit. 



dt ^ 



N 



TV, Uir,)+J2v{\r,^v,\) +r,.(t), 



(2) 

with the pair interaction potential between the particles V{r) 
and an external potential U{r). Both the external potential 
C/(r) and the flow field u(r) can depend on time. For clarity 
of notation we will not make this dependence explicit in the 
equations. The flow field is not necessarily divergence free 
(i.e., the solvent can be compressible). V; denotes the gra- 
dient with respect to r^. We approximate the noise generated 
by the thermal motion of the solvent particles by a Wiener 
process. 



{mm = 



and 



(3) 
(4) 



with the temperature T measured in units of energy (setting 
fee = 1) and the mobility coefficient F > 0. The mobility co- 
efficient has to appear in the correlation in Eq. (|4|i in order to 
fulfill the fluctuation-dissipation theorem and to get the cor- 
rect equilibrium distribution for u(r) = 0. The boundaries of 
V are impermeable for the particles or periodic (or a mixture 
of both), and therefore the number of particles is conserved. 

The Fokker-Planck equation corresponding to the Langevin 
equation (O gives the time evolution of the probability density 
W{ri , . . . ,rN,t) for finding the particles at time t at the po- 



sitions ri, 

dW _ 
~dt ~ 




V,C/(rO 



W 



(5) 



For a potential flow, the velocity field can be written as the 
gradient of a scalar field, u(r) = — F V$(r) ifTsIl, such that 
the external potential and the effect of the flow field can be 
combined into a modified external potential U* (r) —U{y) + 
^(r). If U* is time-independent, one can find a stationary 
probability density Vl^q*(ri, . . . , rjv) which fulfills the detailed 
balance condition for Eq. dU, i.e., the term in curly brackets 
is zero for each i = 1, . . . , iV, 



N 



V,;t/*(r,) + ^ V,V{\r, - r, I) + T V, 



W^q* = 0. (6) 



The solution is 



C/*(r.)+E V(\v,-T,\) 



normalized with the sum of states Z* such that 



d^ri...SrNW*{vi,. 



1. 



(7) 



(8) 



For such a situation the whole apparatus of equilibrium statis- 
tical mechanics can be used in order to calculate expectation 
values and correlations in a stationary non-equilibrium situa- 
tion. But this is restricted to cases when the detailed balance 
condition holds, which implies necessarily a potential flow. 
Even in the Stokes flow ([T]i or in a simple shear flow (e.g., in 
Couette or Poiseuille flow) this is not true because V x u 7^ 0. 
For flows with a finite vorticity there is no detailed balance in 
a strict sense, see [jgl Eq. (5.3.4(c))]. 

From Eq. (|5]l we can calculate the time evolution of the 
noise averaged particle density p{Y, t), namely, the expecta- 
tion value of the density operator /5(r, i) = X]i=i '^('^~'^«(^))' 



| + v.(pu)^v.r 



pVU + TVp (9) 

d3r'y(|r-r'|)p(2)(r,r',t) , 



with the mean density 

p(r, t)^N J J d\2 ■ ■ ■ d?rN W{y, Y2,...,yn, t), (10) 

yjv-l 

and the non-equilibrium density-density coiTelation function 
X yy"dV3...dVwiy(r,r',r3,...,rjv,t). (11) 

yJV-2 
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Eq. (|9]l is the starting point of a hierarchy of N evolution equa- 
tions which connect the time derivative of the n-point density 
correlation function to the n+l-point density correlation func- 
tion, similar to the BBGKY hierarchy for deterministic sys- 
tems with inertia or the BGY hierarchy for equilibrium corre- 
lation functions. 

In order to find a closed equation for the time evolution 
of p{r, t) we approximate the interaction term in Eq. (|9]l by 
its value in an equilibrium system with the same interac- 
tion potential V{r). Let us first restrict our considerations 
to the case that detailed balance holds (so that, in particular, 
u = — r V<i>). We modify our system of Brownian particles 
by applying an external potential ^' (r) so as to create a system 
whose equilibrium density distribution is /5^(r) = p{^,t)- 
This new potential 5'(r) depends on p{r, t) and will be differ- 
ent for each t as long as p{r, t) is not stationary. The Fokker- 
Planck equation for the modified system will be Eq. (|5]l but 
with U{r) replaced by U{r) + ^'(r). The equilibrium proba- 
bility density M^* of the modified system is given by Eq. Q 
but with U*{r) replaced by U*{r) + ^'(r). If we integrate the 
detailed balance condition (|6j for M^* over N — 1 positions, 
we get 



Peq- ^ill '^l- C3 a closure to the hierarchy of 
equations starting with Eq. ©I Hereby we assume that the 
density correlations at time t in the non-equilibrium system 
with mean density p{r, t) are the same as in an equilibrium 
system with the additional potential and with equilibrium 
mean density yO^(r) = p(r, t). We then get 



dp 

at 



V-(pu) 



(16) 



with the free energy functional 



^[p]=^x[p] + y"dV {Tp[ln(pA^) -1] +pC/}. (17) 



In thermodynamic equilibrium for u = and time- 
independent U, the equilibrium density distribution given by 

is a stationary solution of Eq. ( fT6] l. An H- 



^ ^ Sp 

theorem 



dt 



d^rT[p] 



dPrV p 



5p ) 



<0 (18) 



UPeq = r 



Pe*v((7 + *) + rvp* 



+ V / d\'V{\v-v'\)pf{v,v') 



(12) 



("2") 

with the equilibrium pair correlation function p\, {y,y') for 
the modified system in the external potential 5". From equi- 
librium density functional theory one knows that the equilib- 
rium density distribution in the grand canonical ensemble is 
the minimum of the grand canonical functional 



+ j d\ {Tp [In (p A3) - 1] + (C/ext - p) p] 



(13) 



with the thermal wavelength A, the chemical potential /i and 
the sum of all external potentials C/ext = C/ + $ + 5*. The 
excess free energy J^^ [p] summarizes the effect of the particle 
interactions and it is not known exactly in general. We take 
the gradient of the Euler-Lagrange equation following from 
the functional ( fT3] i: since in thermal equilibrium the chemical 
potential is constant across the whole system we get 



[p] 



Sp 



If we compare Eq. (fT2l i with (fT4l) we can see that 

SJhx[p] 



V J dVy(|r-r'|)p(,^)(r,r') = Pe*V 



Sp 



(14) 



(15) 



Note that the right hand side does not depend on the velocity 
potential $ while the dependence on enters only through 



guarantees that the time evolution actually converges to the 
equilibrium distribution. (The dynamics in Eq. (O together 
with the boundary conditions taken for V imply that the par- 
ticle current through the system boundaries is zero and there- 
fore surface terms from partial integration vanish.) The final 
chemical potential is then determined by the conserved num- 
ber of particles in the system. As discussed above, a system 
in a potential flow corresponds to an equilibrium system with 
a modified external potential U* = U + . Eq. ( fT6l ) can then 
be written in the form 



dp 

at 



V TpV 



SJ'*[p] 
Sp 



(19) 



with the modified free energy functional J-*[p] = J-[p] + 
J d^r $ p. Thus we have an H-theorem for T* [p] instead of 

V 

!F[p] and the "equilibrium state" is determined by '^'^^^^ = 

p*. 

The right hand side of Eq. (fTSl l is completely independent 
of the flow field and one could be tempted to use it as a clo- 
sure to Eq. (|9]l for the most general case (i.e., non-potential 
flows). However, this would mean approximating density cor- 
relations in a driven non-equilibrium system where detailed 
balance cannot be achieved by thermal equilibrium correla- 
tions. While we could argue that close to equilibrium Eq. (fTSl l 
may be a reasonable approximation if detailed balance still 
holds, there is no such argument for the most general case that 
detailed balance is violated. The study of Ref. addresses 
a situation where the approximation is found to be good al- 
though there is no detailed balance since a net particle current 
is driven through the system. In the next Section we consider 
some examples with the purpose of assessing (i) the effect of 
a more realistic Stokes flow as discussed in the Introduction, 
and (ii) the validity of the approximate Eq. ( fT6] l for a non- 
potential flow. 
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III. EXAMPLES 

As an example we study a solution of polymers (radius d, 
density po) in an incompressible Newtonian solvent flowing 
around a spherical colloidal particle (radius R) in a stationary 
situation. We model the polymer coils as point-like particles 
from the point of view of the solvent, but with a finite interac- 
tion range cr concerning other polymer coils and D — d + R 
concerning the colloidal particle, the interaction with the lat- 
ter being of hard-wall type, see Fig. [T] The velocity field of 
the solvent is given by the Stokes flow, Eq. ([T]i, and we choose 
c = ce^. Measuring lengths in terms of D, the dimension- 
less parameters determining the system are the Peclet number 
c* = cD/{TT), the colloid radius R* = R/D, and the poly- 
mer's mutual interaction range a/D. 

A. Ideal polymers 

For ideal solute particles in an incompressible solvent with 
U = 0, the stationary condition §f = from Eq. ( fTSI l reads 

u-Vp = TTAp, (20) 

where F T is the diffusion constant of the solute particles. 
The hard interaction with the colloid is written as a bound- 
ary condition for the current density of the solute particles 
j = up-FTVpatr/D = 1, 

{er-i)l/D^,=0. (21) 

We expand the density field p(r) in spherical harmonics up to 
order N and obtain a system of + 1 ordinary differential 
equations for the |r|-dependent coefficients which we solve 
numerically with AUTO 2000 AUTO 2000 is a software 
which solves autonomous boundary value problems for sys- 
tems of ordinary differential equations by continuation, i.e., 
by starting from a known solution for a specific set of prob- 
lem parameters (for c = we have p ^ po) and changing 
parameters (in our case c) continuously until the desired value 
is reached. 

As demonstrated for i?* = and R* =0.9 in Figs.|2ja) and 
(b), respectively, the bow wave effect is large when the colloid 
is small compared to the polymers and does not distort too 
much the flow, reaching its maximum as R* 0. This is the 
case investigated in f^. When the colloid is large compared 
to the polymers, the bow wave effect is small. In the limit 
d/R — > 0, the effect vanishes completely since the polymers 
behave like solvent molecules. Figs. [Sja) and (b) shows the 
density right in front of the colloid as a function of R* and 
of c*, respectively. The density of ideal solute particles scales 
almost linearly with the velocity c* . 

B. Gaussian polymers 

Here we address the case of interacting polymers. We con- 
sider the same polymer-polymer interaction potential studied 




z/D zlD 



FIG. 2: Contour plots of the density of ideal polymers for a flow 
velocity c* = 10. The white circle at the origin is the colloidal 
particle with radius R, the black circle is the annulus of thickness d 
and outer diameter D which is unaccessible to the polymer centers 
due to the hard-wall interaction, c.f. Fig. [T] (a) corresponds to a 
uniform flow u(r) = —ce^ (i.e., to R* — 0). The maximum density 
in front of the colloid is p{r)/po = 6.31. (6) corresponds to R* — 
0.9. The bow effect is reduced drastically. The maximum density in 
front of the colloid is p(r) / po = 1.05. 
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FIG. 3: Density of ideal polymers at the point x = y — 0, z/D — 1, 
i.e., right in front of the forbidden zone around the colloid, (a) shows 
p(0, 0, D) as a function of the colloid size for c* = 1 and c* = 10. 
(b) shows p(0, 0, D) as a function of c* for different values of R* . 

in Ref. lH, namely 

y(r) = rexph(r/a)2]. (22) 

The interaction of the polymers with the colloidal particle is 
modelled as an external potential of the form 

C/(r) = 10rexp[-(r/a)^]. (23) 

This potential rises steeply up to 10 T, thus resembling a hard 
wall. The length a must be related to the radius of the for- 
bidden zone D around the colloidal particle. We convention- 
ally set the value of a by the condition U{D) = 2T, giv- 
ing a « 0.924 Z?. We take a = 2d (i.e., additive mixture of 
polymers and colloidal particle) and R = 1.7 a, leading to 
R* w Q.77, (j/D « 0.46. Finally, we also considered the 
choice R* = 0, cr/D w 0.46, which represents a hard particle 
that does not distort the uniform flow (i? = in Eq. ([T]i) in a 
non-additive mixture (a ^ 2d): this was the model addressed 
in Ref. fl. 

We ran Brownian dynamics (BD) simulations of this system 
for two values of the flow velocity corresponding to the poly- 
mer Peclet numbers {a/D)c* = 1 and 10 studied in Ref. yll 
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(i.e., c* PS 2.2 and 22). We considered a colloidal parti- 
cle at the center of a box of dimensions ~ Ly = 12 a 
and Lz = 24 tr with periodic boundary conditions. The box 
contained N = 3456 polymers, corresponding to a mean 
polymer number density pocr^ = 1. We took a timestep of 
0.003 (T^r/T for the discretized Langevin dynamics. The sys- 
tem was allowed to relax for 10^ timesteps, after which collec- 
tion of data was carried out during 10^ timesteps. Even though 
the simulated system is finite we used the analytically known 
flow field around a sphere in an infinite medium, Eq. ([TJ. The 
error due to the truncation of this flow by the boundary of the 
simulation box is largest (about 20%) at the midplane of the 
colloid [z = 0). This introduces effectively a discontinuity in 
the flow velocity field at the boundary which we discuss later 
We also solved numerically the dDFT in the random phase 
approximation (a mean-field model), i.e., with 

^x[p] = ^JJ d'rd'r' Vi\r - r'|)p(r)p(r'). (24) 

The time evolution given by Eq. ( fT6b of an initially homoge- 
neous density was solved in cylindrical coordinates on a grid 
spanning the domain —120 < z/a < 24, < r±/(j < 60, 
where r± = y^x^ + y^. The grid constant was 0.0125 a near 
the colloid, i.e., for rj^ < 6 cr, and 0.1 cr in the rest of the 
domain. For details on the numerical procedure see The 
boundary condition at the domain border was p ~ po also in 
this case we used the flow field given by Eq. The error in- 
troduced here is smaller than in the BD simulations since the 
integration domain for the dDFT is larger than the BD simu- 
lation box. 

Fig. m presents the density field p{z), spatially averaged 
over thin disks of radius a and thickness 2 Az = 0.05 a cen- 
tered at the z-axis, i.e., 

z-\-Az cr 

p{^) ■= ^2^^ J dz' J dr±r±p{r±,z'). (25) 

z-Az 

The results of both the BD simulations and the dDFT illus- 
trate the dramatic effect of advection by the Stokes flow ([T]). 
In particular, at the higher velocity c* = 22 and R* = 
(uniform flow) there is a marked accumulation of polymers 
in front of the particle and a strong depletion behind it which 
are hardly observable for R* = 0.77. In general, the effect 
of the Stokes flow is to weaken the influence of the colloid 
on the density profile, as the polymers tend to be advected 
by the stream and to travel around the particle. Actually, for 
c* = 2.2 (and smaller) the deformation by the Stokes flow is 
so tiny that the dDFT profile p{z) in Fig. Ufa) is indistinguish- 
able from the equilibrium profile (i.e., c* = 0). 

We notice a discrepancy between the density profiles mea- 
sured in the BD simulation and those calculated numerically 
in the dDFT. We attribute this to two finite-size effects in 
the simulation which have been confirmed by performing BD 
simulations in a smaller box at the same polymer number den- 
sity (L^ = Ly = 8a, Lz ^ 16c7, N = 1024). First, 
the numerical solution of the dDFT for a Stokes flow with 
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FIG. 4: Plots of p{z) defined in Eq. i25\ as provided by the nu- 
merical solution of dDFT (lines) and as measured in BD simulations 
(symbols), (a) corresponds to a velocity c* = 2.2 and (b) to c* = 22. 
In each plot the results for both uniform flow (R* — 0) and Stokes 
flow (R* = 0.77) are presented. 



c* = 22 exhibits a long (« 15 a) tail of slight polymer de- 
pletion (p « 0.98 po) behind the colloidal particle. The tail is 
longer than the length of the BD simulation box in z-direction 
and, as a consequence of the periodic boundary conditions, 
the inflowing density far ahead of the particle is smaller than 
Po- This screening effect is very noticeable when the flow is 
approximated as uniform because the depletion of polymers 
behind the colloidal particle is very large, see Fig. |Ub). But 
in the case of the Stokes flow, this effect seems to less impor- 
tant compared to the second effect: the discontinuity of the 
normal component of the flow field at the lateral boundaries 
of the simulation box leads to a non-vanishing divergence of 
the flow there (we remind that the flow described by Eq. ^ 
has V • u = 0). Fig.|5]represents the density profile averaged 
over thin cylindrical shells of height a and radial thickness 
Ar = 0.05 a coaxial with the z-axis, i.e., 

P{r±,Zc) , ArN A / / <^^'±T'±Pi^'±^^')- 

[r± + —) Ar J J 

Za-§ r± 

(26) 

The periodic boundary conditions imply that V ■ u < effec- 
tively at the side boundaries located upstream, where there- 
fore the density is enhanced: Even though we expect the den- 
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FIG. 5: The plot represents /5(r^,2c), see Eq. ( I26t . in a Stokes 
flow (IV — 0.77) with c* = 22 upstream (at Zc = 1.5a, triangles 
up) and downstream (at Zc = —1.5a triangles down) of the colloid, 
as measured in BD simulations. The inset shows the results for a 
smaller simulation box. 

sity to decrease towards po as the radial distance r± to the 
colloid increases we find instead an increase of the density at 
the boundary of the simulation box. At the side boundaries lo- 
cated downstream, on the other hand, V • u > and the region 
near the boundary of the simulation box becomes depleted of 
polymers. As expected this effect is enhanced for reduced 
box size, see the inset in Fig. |5] The comparison of the BD 
results with the dDFT calculation in Fig. |4]indicates that the 
overall consequence of these effects is a density enhancement 
near the colloidal particle. In view of these important finite- 
size effects, we cannot quantify the validity of the approxima- 
tion given by Eq. ( fTSI l for a realistic flow. When there is only 
uniform flow [R* = 0), however, the finite-size effects are 
much less pronounced and we find a good agreement between 
BD simulations and dDFT calculations, in concordance with 
Ref. il. 



C. Drag force on colloids 

We have also measured the force in the z-direction exerted 
by the polymers on the particle. This force is additional to the 
Stokes drag force Fstokes exerted by the flowing solvent. If F in 
Eq. (|2]i is assumed to be given by the Stokes-Einstein relation 
for spherical polymers of diameter cr, the Stokes drag for a 
colloid of radius R in the same solvent is given by Fstokes = 
2 c* {R/cr) (T/D). For example, for the colloid radius R* = 
0.77 considered in the BD simulations this gives Fstokes = 
3.4 c* T/D. 

In the BD simulations, the force excerted on the colloid 
by the polymers can be measured directly. In the case of 
ideal particles discussed in Sec. lIII Al we use the ideal gas law 
p = T p in order to calculate the local pressure on the colloid 
surface. Integrating the local pressure over the surface yields 
the force on the colloid. 



Type of flow 


c* 




ideal 


-fstokcs 


uniform 


2.2 


42.5 


41.9 




Stokes 


2.2 


6.38 


3.14 


7.48 


uniform 


22 


215 


290 




Stokes 


22 


10.3 


18.0 


74.8 



TABLE I: Mean force exerted by the polymers measured in the BD 
simulations for different types of flow (R* = for the uniform flow, 
and R* = 0.77 for the Stokes flow), compared to the force exerted 
by ideal polymers and to the Stokes friction of the colloid. The forces 
are given in units of T/D. 



simulation box size 


c* 






small 


2.2 


11 


63.4 


large 


2.2 


6.4 


103 


small 


22 


19 


65.6 


large 


22 


10 


103 



TABLE II: Mean and variance of the force exerted by the polymers 
in a Stokes flow measured in BD simulations of different boxsizes 
(see Sec. nil BV The forces are given in units of T/D. 

Table H] collects the mean force for different types of flow 
(R* = for uniform flow and R* — 0.77 for Stokes flow) 
and velocities. The results confirm the necessity to take the 
solvent flow into account: the mean force in the case of Stokes 
flow is markedly smaller (at most of the order of Fstokes) than 
in the case of a homogeneous flow and the dependence on c* 
is milder. This can be understood in terms of the reduction of 
the bow wave effect in the density profile around the particle 
by the Stokes flow in the solvent which advects the polymers. 
The forces in the BD simulations are of the same order of 
magnitude as in the ideal case, but in the case of the Stokes 
flow they have a weaker dependence on c* . 

Because the sizes involved are of the order of the micro- 
scopic length D, the variance of the force measured in the 
Brownian dynamics is relatively large. However, we find that 
it is not affected by the flow type and velocity and it coincides 
with the variance of Fz we have measured in the equilibrium 
state (c* = 0). For comparison, Tab.HJcoIIects the force mea- 
sured in the BD simulation in the smaller box. The increased 
force is consistent with the density enhancement near the par- 
ticle caused by the finite-size effects. We also attribute the 
increase of the fluctuations to these effects. 



IV. CONCLUSIONS 

We have proposed a dynamic density functional theory 
(dDFT), Eq. ( fT6l ), for interacting Brownian particles in a flow- 
ing solvent under the assumption that detailed balance holds 
(which requires, in particular, a curl-free flow). We get the 
same equation as already derived in ^ but with the partial 
time derivative replaced by the total (material) time deriva- 
tive. The whole effect of the flow field can be summarized 
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into a modified external potential, allowing application of the 
whole machinery of equilibrium statistical mechanics. Thus, 
we are able to find an H-theorem for a modified free energy. 

In this paper we include the displacement of the solvent 
by the colloid, but the hydrodynamic interactions between the 
colloid and the solute particles as well as the hydrodynamic 
interactions among the solute particles were not taken into ac- 
count. While the latter is a highly non-trivial and still open 
problem, the first can be included in a straightforward manner 
by replacing the mobility F by a space dependent and sym- 
metric mobility tensor r(r,i) in Eqs. (|2]i and (|4|i. Thereby the 
noise becomes multiplicative and the appropriate calculus has 
to be considered such that it leads to the Fokker-Planck equa- 
tion (O with r replaced by T{ri). Then the equilibrium dis- 
tribution is not changed. For spherical particles in the vicinity 
of planar walls the mobility tensor can be calculated in the 
limit of large distances ifTTll . This result has been extended to 
surfaces with a partial slip boundary condition in lfl2ll . Due 
to the translational symmetry of the system, F is diagonal. 
While the mobility perpendicular to the wall increases with 
distance, the distance dependence of the mobility parallel to 
the wall depends on the slip condition. For no-slip it increases 
while for total slip it decreases with the distance to the wall. 
The hydrodynamic interaction between two spheres has been 
calculated, e.g., in the Rotne-Prager approximation ifisi [l4ll . 
The derivation of Eq. ( fT6] l essentially carries through with the 
only exception that the quotient u(ri)/r in Eq. (|5]l has to be 
replaced by a field u(rj) with T{ri) ■ u(ri) = u(ri). In or- 
der to absorb the flow field into a modified external potential, 
u (and not only u(ri)) has to be curl free with u = — V$. 
Instead of Eq. ( fTSI l we then get 



In Ref. ill], the polymer distribution was studied in a poly- 
mer solution flowing uniformly through a spherical particle 
which is hard only for the polymers. In spite of the violation 
of detailed balance by the boundary conditions, the compar- 
ison between simulations and the numerical solution of the 
proposed dDFT was good. In this paper we have considered 
the more realistic case of a Stokes flow ([T]i around the parti- 
cle. The exact solution of the ideal case (no polymer-polymer 
interaction), the numerical solution of the interacting case as 
well as the corresponding Brownian dynamics simulations ev- 
idence all the dramatic effect by advection on the properties of 
the stationary solution. We conclude that the approximation of 
uniform flow, as employed in Refs. iH 0,111], is quantitatively 
bad. We have found discrepancies in the density distribution 
of polymers as measured in the simulations and as computed 
numerically in the framework of the dDFT. However, the dis- 
crepancies could be rationalized in terms of finite-size effects 
in the simulations due to the slow decay of the Stokes flow far 
from the obstacle. Thus, although a quantitative check of the 
validity of the approximations leading to the dDFT in Eq. ( fTSI l 
and of its validity for non-potential flows was not possible the 
results are encouraging. 
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